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Abstract 

In this paper we are interested with a strongly coupled system of partial differential 
equations that modelizes free convection in a two-dimensional bounded domain filled with 
a fluid saturated porous medium. This model is inspired by the one of free convection 
near a semi-infinite impermeable vertical flat plate embedded in a fluid saturated porous 
medium. We establish the existence and uniqueness of the solution for small data in some 
unusual spaces. 



1 Introduction 

In the literature, many papers about free convection in fluid saturated porous media study the 
case of the semi-infinite vertical flat plate in the framework of boundary layer approximations. 
This approach allows to introduce similarity variables to reduce the whole system of partial 
differential equations into one single ordinary differential equation of the third order with ap- 
propriate boundary values. This two points boundary value problem can be studied using a 
shooting method or an auxiliary dynamical system either in the case of prescribed temperature 
or in the case of prescribed heat flux along the plate. 

In this article we first present the derivation of the equations, show how the boundary layer 
approximation leads to the two points boundary value problem and the similarity solutions, 
then we rewrite the full problem of free convection in a two-dimensional bounded domain filled 
with a fluid saturated porous medium. This new model, written in terms of stream function 
and temperature, consists in two strongly coupled partial differential equations. We establish 
the existence and uniqueness of its solution for small data. 



AMS 2000 Subject Classification: 34B15, 34B40, 35Q35, 76R10, 76S05. 
Key words and phrases: Free convection, porous medium, coupled pdes. 
f s.akesbi@uha.fr I b.brighi@uha.fr tj j-d.hoernel@wanadoo.fr 



1 



2 The semi-infinite vertical flat plate case 



Let us consider a semi-infinite vertical permeable or impermeable flat plate embedded in a 
fluid saturated porous medium at the ambient temperature Too, and a rectangular Cartesian 
co-ordinates system with the origin fixed at the leading edge of the vertical plate, the x-axis 
directed upward along the plate and the y-axis normal to it. If we suppose that the porous 
medium is homogeneous and isotropic, that all the properties of the fluid and the porous 
medium are constants and that the fluid is incompressible and follows the Darcy-Boussinesq 
law we obtain the following governing equations 

du dv q 
dx dy 

k f dp 

u = — 77- + pg 



p \dx 

k dp 

v = 7T' 

pdy 

dT dT , (d 2 T d 2 T 
a— +v— = A — r + 



dx dy \ dx 2 dy 2 

p = p 00 (l-P(T-T 00 )) 

in which u and v are the Darcy velocities in the x and y directions, p, p and (3 are the density, 
viscosity and thermal expansion coefficient of the fluid, k is the permeability of the saturated 
porous medium, A is its thermal diffusivity, p is the pressure, T the temperature and g the 
acceleration of the gravity. The subscript oo is used for values taken far from the plate. In 
our system of co-ordinates there are two main interesting sets of boundary conditions along the 
plate. 

First, the temperature is prescribed on the wall that gives 

v(x,0) =ujx m ^ 1 , T(x,0) = T w (x) =T oa + Ax m (2.1) 

with m e R and A > 0, see HE], [IB!, EI], [2H] and [S2]. 

Secondly, the heat flux is prescribed along the plate that leads to 

v(x, 0) = wi't 1 , — (x, 0) = -x m (2.2) 
dy 

with m G R, see PS] and [TTj . 

The parameter u G R is the mass transfer coefficient. For an impermeable wall we have 
u = 0, and for a permeable wall, u < corresponds to fluid suction and uj > to fluid injection. 
The boundary conditions far from the plate are the same in both cases (jH.Hj) and (|H.4jl 

u(x, oo) = 0, T(x, oo) = T^. 

If we introduce the stream function ^ such that 

_ d^ <9# 

dy ' dx 
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we obtain the system in which it remains only \& and T 



c> 2 * | _ Poo /3gk dT 
dx 2 dy 2 // dy ' 

A (cPT_ + _ dTd^_ _ dTd^_ 

\ dx 2 dy 2 J dx dy dy dx 

Along the wall, the boundary conditions ()3.3|) become 

Oty m-l 

— (x,0) = -ux— , T(x,0) =T w (x) =T OQ + Ax m (2.5) 
ax 

and (|3.4|) becomes 

m-l <9T 

— (x,0) = -o;x-i- I — (*,0) = -*"\ (2.6) 
The boundary conditions far from the plate become 

— (x, oo) = 0, T(x, oo) = T^. (2.7) 

We will start from the equations ()2.3j) - (j2.4j) subjected to the boundary conditions (|2.5jl and 
(12. 7|) with a; = to write a new model, settled in a two-dimensional bounded domain, that we 
will study in the rest of this paper. 

Before doing this, let us say a few words about the similarity solutions. Assuming that con- 
vection takes place in a thin layer around the plate, we obtain the boundary layer approximation 

d 2 ^ _ Poo (3gk dT 
dy 2 fi dy 

d 2 T _ 1 (dTd^ dTd^\ 
dy 2 X \dx dy dy dx J 

with the same boundary conditions (|2.5jl or (|2.fij) and (|2.7jl as before. 

For the case of prescribed heat, introducing the new dimensionless similarity variables 

t=(Ra x )^, *(x,y) = \(Ra x )*f(t), T(x,y) = (T w (x) -T^t) + T^ 
x 

with 

p ao Pgk^T w {x)-T 00 )x 

Ka x — 

/iA 

the local Rayleigh number, equations ()2.8|) and ()2.9|) with the boundary conditions ()2.5|) and 
()2.7|) leads to the third order ordinary differential equations 

f" + ^ff" - mf = 

on [0, oo) subjected to 

/(0) = -7, /'(0) = 1 and /'(oo) = 



(2.9) 
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where 

7 = ^ A / I ", 
m + 1 y PooflgkAX 

One can find explicit solutions of this problem for some particular values of 7 or m in jj], 
0, [0], |20], pp] - [2"H] . PO] and [33]. For mathematical results about existence, nonexistence, 
uniqueness, nonuniqueness and asymptotic behavior, see |2j, [E], [Hj and |2H| for 7 = 0, and [H], 
[12], [U3, [23] and j21] for the general case. Numerical investigations can be found in [2], [Zj, 
HE], HE], [2B], |SH] and [SB]. 

In the case of prescribed heat flux, we introduce the new dimensionless similarity variables 

t = 3-^Rlx n ^ 1 y, *(x,y) = 3%Rhx m ¥f(t), T{x, y) = £^x^6{t) + 
and the Rayleigh number 

H-a — : — • 

Then, equations ()2.8|) and ()2.9|) with the boundary conditions (I2.6j) - (j2.7|) give 

r + (m + 2)/r-(2m+l)f 2 = 

and 

/(0) = -7, /"(0) = -l and /'(oo) = 

where 

1 —i 
7 ~ A(m + 2) ' 

The study of existence, uniqueness and qualitative properties of the solutions of this problem 
is made in [TO] . For a survey of the two cases, see [H]. This equation can also be found in 
industrial processes such as boundary layer flow adjacent to stretching walls (see [2], |S], (20], 
[20], |30j) or excitation of liquid metals in a high-frequency magnetic field (see [33]). 

One particular case of the two previous equations is the Blasius equation /'" + //" = 
introduced in [S] and studied, for example, in [1], [TO] and [2*7] . 

The case of mixed convection /"' + //" + mf'(l — /') =0 with m G R is interesting 
too and results about it can be found in pQ, ^3], [2E] and [31] • The Falkner-Skan equation 
/'" + //" + m{l - f' 2 ) = with m E R is in the same family of problems, see [U], [22], [23, 
[20], EH], [33 and [10] for results about it. 

New results about the more general equation /'" + //" + g(f') = for some given function 
g can be found in [H], see also 

3 A model problem in a bounded domain 

Let Q C M 2 be a simply connected, bounded lipschitz domain whose boundary T = dfl is 
divided in two connected parts Ti and T 2 such that 

fi u r 2 = r and Ti n r 2 = 0. 
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We start from the previous equations ()2.3)) - (j2.4j) in terms of the stream function \1/ and the 

temperature T with K = ^0, poo ^ gk j , and assuming that I\ is impermeable and that the 

temperature T w > is known on the whole boundary T, we modify the equation ()2.3j) by 
setting K(x) = (ki(x), k 2 (x)) G M 2 with < ||i^||oo < oo. Then, we obtain the following new 
problem in the bounded domain Q, which consists in finding (\&,T) 



* : Q 

T : Q, 

verifying the equations in Q 



A\l> = K.VT, (3.1) 
AAT = VT.(V^) ± , (3.2) 



the boundary conditions on T for \1/ 



<9^ 

f = on Ti and — = on T 2 (3.3) 
on 

and the boundary conditions on T for T 

T = T W onT (3.4) 
where A G M + * and for all x = (u, v) G f2, let x 1 - = (v, —u). 

3.1 Preliminary results 

Let us assume that T w G H^(V) and let be the unique function in H 1 ^) verifying 

A0 = in tt, (3.5) 
9 = T w on T. (3.6) 

In the following we will need that VO G L°°(Q), thus we will suppose that it holds (it is the 
case if T w G H2(T) for example). 
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If T) is a solution of ()3.1|) - (j3.4|) and if we set H = T — 6, then H) is a solution of 

Afy = K.VH + K.VQ, (3.7) 
XAH = V#.(W) X + Ve.(V^) 1 " (3.8) 

in the domain f2 with the boundary conditions for \1/ 

f = on and — = on T 2 (3.9) 
on 

and the boundary conditions for H 

F = on T. (3.10) 

Conversly, it is clear that if if) is a solution of (l3~7|l - (j3~Tn|) then T) := if + 6) is a 
solution of (JSZEJ-flOl). 

In the following we set ||.||n(n) = ll-lli, ||-IU 2 (n) = IMk ||-IU°°(n) = IMU and 

(u, v ) = / uvdx. 
Jn 

Definition 1 For u G f G #o(^) and w G if 1 ^) /et 

a(u,v,w) = (iiVj),(V«j) 1 )i2 (!2)ii 2 (!]) . 

Remark 1 TTie trilinear form a is well defined because for u G L°°(Q), v G Ho(Q) and w G 
i/ 1 (fi) we /iat>e 

\a(u,v,w)\ < ||m||oo|| Vv|| 2 || Vw|| 2 - 

Proposition 1 For u G Hq(Q) fl L°°(fi) and u G iJ 1 (fi) we /icwe 

a(n, it, = 0. (3-11) 

Proof. First, let us notice that if u G H^(Q) n L°°(fi) then u 2 G flj(fi) and V(n 2 ) = 2«Vw. 
Hence 

a(u,u,v) = (uVu, (Vt)) i ) L 2 {! j )iL 2 (S]) 
= 2^ Vw2 ' ( Vw ) ± )i 2 (n),i 2 (n) 



1 (div((Vt») ± ), M 2 ) H -i (n)i ^i (c) 
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because u = on T and div((Vf) _L ) = in if -1 . ■ 
Remark 2 Foru,v G i/o(Q) fl L°°(fi) and w G if^fi) we /lave 

a(u, v, w) = —a(v, u, w). (3-12) 
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3.2 A priori estimates 

Let 

W 9 = {u\ue H\Q) and u = on r±} 

and 

w H = H^(n)nL°°(n). 

The spaces Wq? and Wh are equipped with the norms \\-\\w 9 an d defined by 

\\ u \\wv = ||Vm|| 2 and = + \\Vu\\l. 

In the following we will use the notation C for the Poincare's constant of Q. 

Definition 2 We will call {^,H) G Wq, x W H a weak solution of the problem (pOty- pHT^) ?/ 

and only if we have 

(W, Vu) + (K.VH, u) + (K.V6, u) = 0, (3.13) 
A(V#, Vu) + a(u, *) + a(u, 6, *) = (3.14) 

for all u G W<$ and v G Wh ■ 

Proposition 2 Let {^,H) e W^ x W H be a solution of the problem (pH^-PT^l) and T = 



H + Q, then 

inf T w < T < supT„,. (3.15) 



r 



r 



Proof. Set I = sup r T w and T+ = sup(T -1,0). As T+ G W H , using (jHUD with u = T+ and 
noticing that (V6, VT + ) = because A© = 0, leads to 

A(VT, VT+) + a(T + , T, *) = 0. 

Using the facts that A(VT,VT+) = A(VT+,VT+) and a(r+ T,tf) = a(T+,T+,^) = by 
proposition [T] we obtain that 

||VT + || 2 = 

and as T + G -f^o (fi) we have T + = on Q. We proceed in the same way with I' = infr T w and 
T~ = inf (T — 0) for the other inequality. ■ 



Proposition 3 Let {^,H) G Wy x W H be a solution of the problem $3.1fy - ffiJ$ , then for 

e have 

|W|| 2 < 2C||Ar|| 00 ||Ve|| 2 and \\VH\\ 2 < IIVOIU. 



lV0||oo < 2C , 2 || X |i we have 
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Proof. Taking u = ^ in (|3.13|) and using Poincare's inequality we obtain 

|| || | < | (i^. vu, I + 1 . ve, I 
<||x||oo(l|Vff|| a + ||ve|| 2 ) ||*|| 2 

^CWKlUiWVHh + WV&h) ||W|| 2 

and 

|| W|| 2 < (711^^(1^^112+ || V6|| a ). (3.16) 
Taking v = H in ()3.14j) leads to 

X(VH, VH) + a(H, H, *) + a(H, 6, *) = 0. 

Then, by proposition ^ we have 

\\\VH\\l<\a(H,Q,*)\ 

< ||ve|U|#|| 2 ||w|| 2 

< ciiveiuiiv^iuiiv^iu 



using Poincare's inequality and 



||Viy|| 2 <^||ve|unv*|| 2 . (3.17) 



Then, combining (jH.lfij) and ()3.17|1 leads to 



||w|| 2 < cii^iunveih + ^M^||ve|U||v*|| 2 . 



Thus 

c 2 \\k\\ 



1 - — ^^liveiu ) ||v*|| 2 < c||^|U||ve|| : 



and as ° "f 1 |o ° || Ve||oo < 1/2 we have 

||W|| 2 < 2C7|| J K'|| 00 ||VO|| 2 . 
Using this new inequality in ()3.17|) . we obtain 

l|Vif|| 2 < ||V6|| 2 . 

■ 

Remark 3 As 

||V9|| 2 < (mes^^lVeiU and \\VQ\U < 



2C 2 || J RT|| 00 
we can rewrite the previous result as 

||W|| 2 < ^(mes^ and \\VH\\ 2 < ^^- (mes »)i 



3.3 Main results 

Theorem 3 Let M = sup r T„,. If MCWKW^ < A, then the problem PTTft- fpH^) admits at 
most one weak solution in x Wh- 

Proof. Let (*i,i21) and (^ 2 ,H 2 ) be two solutions of (f3~Tjl - (f3~TTHl . Setting H — H x — H 2 and 
v[r = ^ - ^r 2 we obtain 

(W,Vu) + (K.VH,u) = 0, 
A(VJ2, Vu) + a(v, H h ^1) - a(u, #2, #2) + a(u, 9, = 

for u G W 7 ^ and t> G Wh- Choosing u = ^ and v = H leads to 

(V*,V*) + (OT,*) = ) (3.18) 
A(V/2, V/2) + a(#, F 1; *x) - a(F, # 2 , #2) + ©, *) = 0. (3.19) 

From equation f)3. 18j) we deduce that 

||W || 2 < C\\K\U\VH\\ 2 . (3.20) 

Let us compute 

a(H, H 1} *x) - a(F, # 2 , * 2 ) = -a(J2 2 , /2 1; - a(H u H 2 , * 2 ) 

= a(/Ti, J H 2 ,*i)-a( J ffi,fl 2 ,* 2 ) 
= a(H 1 ,H 2 ,V) 
= a{H,H x ,y). 

Thus, using now equation ()3.19|) we get 

A(Vi2, Vi2) + a(#, #1 + ©, f ) = 

and 

\\\VH\\l < \a{H i H 1 + Q,^)\ 

< 10(71, I 

< llTilUHV^IhUV^Ih 

< M||VJ2|| 2 ||Vf || 2 



with M = suppT^. Therefore 



and using (|3.20|) we have 



|V/2|| 2 < -H|W|| 2 
A 



|vfi||, < MC flh|yg||, 

A 
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Choosing j l|o ° < 1 we obtain ||Vi?||2 = and || \7\& || 2 = 0. This complete the proof. ■ 

In the following Theorem, we prove the existence of a strong solution H) of the problem 
(|3.7|) - ([3.10j) under some hypothesis on the data. To this aim, let us define the spaces 

f du 1 

W* = < u | u e H 2 (tt), u = on r\ and — = on T 2 \ 

and 

w H = H^(n)nH 2 (n). 

These spaces are equipped with the following norms 

IMIiW = ll^ll/w 
\\ v \\w H = \\^ v \\m(n), 

and 

||(M,f)|| L 2 (Q)xL 2 (n) = ||u|| 2 + ||u|| 2 . 

Theorem 4 Let M = sup r T w . For max {C|| VO||oo, M } < and small values of \\K.VQ\\2, 

there exists a unique solution (^,H) of the problem \3. TD - fj. 1 ?. ifljl in the space x Wh- 

Proof. Let us define the operator 

A : WyxW H ^ L 2 (Vl) x L 2 (n) 

such that A(V,H) = , H), A 2 (V , H)) with 

Ai(V,H) = A$ - K.VH, 

A 2 (V, H) = XAH - V/f.(V^) ± - Ve.(W) 1 . 

Let us remark that, using the Sobolev embedding theorem, we have •— > L A (Q) in such a 
way that V//.(V^) ± e L 2 (fi). 

In term of the operator A, the equations ([3.7)1 - ([3.8)1 can be rewritten as 

A(V,H) = (K.VQ,0). 

Notice that H) = (0, 0) is a solution of A(^f, H) = (0, 0) and by the same argument as in 
Theorem El it is the only one. 

Now we want to show that the solution of A(^f, H) = (K.VQ, 0) also exists for small values of 
||K.V6|| 2 . To this end, let us compute the Frechet derivative of A. For e and G G Wh, 
we have 

A{<f), G) - (A0 - K.VG, XAG - VG.iV^) := G) - L(<f>, G) 

= (O,-VG.(V0) ± ) 
= G)\\^ x w H ) 
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because 



||(0, VG.(V0) x )|U 2(n)xL2(n) = ||VG.(V0)iU 2( n) 

< ||VG||£4 ( n)||W|| L 4 ( n) 
<C s 2 ||VG|| Hl(n) ||V0||Hi(n) 

<c 2 s \\^g)\\%^ Wh 

where C s is the Sobolev constant corresponding to the continuity of the embedding H l (Vt) "—>■ 
L 4 (fi). Thus, L defined by L(<f>, G) = (Acp — K.VG, A AG- Ve.fV^)- 1 ) is the Frechet derivative 
of A at the point (0, 0), i.e. 

A'(0, 0).(0, G) = (A0 - K.VG, XAG - Ve.fV^). 

For / and g in L 2 (f2) let us now consider the system A'(0, O).(0, G) = (/, (?) that can be written 

as 

-A</> + K.VG = f, (3.21) 
-XAG + VQ.(V(f)) L = g. (3.22) 

To prove the existence of a solution (\l/, H) of ()3.7j) - (j3.1())) it remains to show that the linear 
operator A'(0, 0) : W 7 ^ x W# — > L 2 (fi) x L 2 (Q) is invertible. To this end, we must first prove 
that for every given / and g in L 2 (Q) the system (j3.21|) - ()3.22|) admits at least a solution and 
secondly that for (/, g) = (0, 0) only (0, G) = (0, 0) is a solution of (f3~2T) - (f3~22]l . 

• First, we want to prove that for every given / and g in L 2 (f2) the system (|3.21|) - (j3.22|) 
admits at least a solution. To this aim, let us define the operator T = Q o S : G i— » Gi 
from i7 1 (Q) into H 1 ^) with S : G ^ <p where <p is the solution of 

-A(f> + fC.VG = / 

in Q with the boundary conditions <p = on T\ and |^ = on T2, and Q : <fi t— > Gi where 
Gi is the solution of 

-AAGi + Ve.(V0) ± = # 
in with the boundary conditions G\ = on T. 

Suppose now that G and G' are given in H X (Q). Let us consider = S(G), (f)' = S(G') 
and G 1 = Q((p), G[ = Q{(j>'). Setting G = G - G', = - </>' and G x = G 1 - G[, by 
(|3.21)l - (|3.22jl we have the inequalities 



\\V(j>\\ 2 dx= - \ (K.VG)4>dx<C\\K\\ 00 [ / \\V(f)\\ 2 dx) ( / \\WG\\ 2 dx 

Jn 

and 

A / \\VG 1 \\ 2 dx = - I VQ.iV^G^x < GUVeHoo ( / ||V0|| 2 ds) ( / {{VG^dx 
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Combining these two inequalities, we obtain 



C 2 ||iTIUIve 

II VG'i||i2(Q) < 

that shows us that if 



VGi|| L 2 ( q) < ^ — ||VG||l2(q) 



x 

then T is a contraction from H l (Q) into itself and admits a fixed point G e H 2 (Q) that 
gives us a solution (0, G)ef$x W/ H of (|S2H)-(|SI221)- 

• The system Ij3.21|) - (|3.22j) with (/, g) = (0,0) admits (0,0) for solution, let us show that 
this solution is unique. Let us suppose that (0, G) G Wq, x Wh is a solution of ()3.21jl - 
(13.22)) . multiplying ()3.21j) by 0, (13.22)1 by G and integrating on Q leads to 

/72||z^|| ||\7ftll 

|)VG|| a < - 11 ll ^ l|VU||oo ||VG'|| 2 
A 

from which we deduce G = and = if C 2 1| 7^ || oo || V© || oo < A. 

This shows that, for small values of ||iT.V©||2, the problem A(ty, H) = (iT.VO, 0) does have 
solutions. Thus, for such values of and K and C 2 || if||oo|| VO||oo < A, the problem (j3.7j) - (j3.10j) 
admits at least one solution (0, G) in Wq, x Wh and, as Wq, x Wh C Wq, x Wh, by Theorem 
EI it is unique if, in addition, we have MCH-KTHoo < A. ■ 

Remark 4 Since, in the previous Theorem we have 

||#.ve||2 < ||K||oo||ve|| 00 (mes n)i 

and 

llveiL < 



C 2 \\K\ 



the condition ||iT.VO||2 small is realized when ^(mes VL)% is small. It is the case, for example, 
when the domain Q is large and the parameter X, that is the thermal diffusivity of the porous 
medium, is small. 

Corollary 1 Let T w e H^(V) and M = sup r T„,. If max {C\\ VO||oo, M} < qw^ti there exists 

a unique solution T) of the problem (Iff. ijl - Qff.^l ) in the space Wq, x H 2 (Q) for small values 
o/||KVe|| 2 . 

Proof. It follows immediately from Theorem 0] and the fact that problems ()3.1j) - ()3.4j) and 
()3.7|) - ()3.10)) are equivalent. ■ 
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4 Conclusion 



In this paper, starting from the model of free convection in a fluid saturated porous medium near 
a semi-infinite vertical flat plate we have written an extension describing this phenomenon in a 
two-dimensional bounded domain. This new problem is given by two strongly coupled partial 
differential equations, that allows us to compute the stream function and the temperature of 
the fluid in the porous medium. 

In a first approach of this complex problem, we have proved existence and uniqueness 
of a solution for small data when a part of the boundary of the domain is assumed to be 
impermeable. 
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